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Abstract 

Inverse problem to recover the skew-self-adjoint Dirac-type system 
from the generalized Weyl matrix function is treated in the paper. 
Sufficient conditions under which the unique solution of the inverse 
problem exists, are formulated in terms of the Weyl function and a 
procedure to solve the inverse problem is given. The case of the gen- 
eralized Weyl functions of the form 0(A) exp{— 2i\D}, where is a 
strictly proper rational matrix function and D = D* > is a diagonal 
matrix, is treated in greater detail. Explicit formulas for the inversion 
of the corresponding semiseparable integral operators and recovery of 
the Dirac-type system are obtained for this case. 
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1 Introduction 

The skew-self-adjoint Dirac-type system 

— u(x, A) = (iXj + jV(x)}u(x, A), x > 0, (1.1) 

where 
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I p is the p x p identity matrix, and v is a p x p matrix function, is actively 
studied in analysis and soliton theory (see, for instance, [H H2] and the ref- 
erences therein). System (11 .11) is an auxiliary linear system for the focusing 
matrix NLSE, sine-Gordon and other important integrable equations. 

The inverse problem to recover a self-adjoint Dirac type system from its 
Weyl or spectral function is closely related to the inversion of the integral 
operators with difference kernels, see [HI [2H1 E21 EES EZj and various references. 
For the discrete analogues of Dirac systems, Toeplitz matrices appear instead 
of the operators with difference kernels [7J [TOj [151 EH] • (Various results on 
Toeplitz matrices and related j-theory one can find, for instance, in [HI [HI dSl 

When the Weyl functions of the self-adjoint Dirac type system are ra- 
tional, one can solve the inverse problem explicitly. One of the approaches 
to solve the inverse problem explicitly is connected with a version of the 
Backlund-Darboux transformation and some notions from system theory 
[221 122] ■ (See also M, M, W\ for this approach, and see |3H] and the 
references therein for explicit formulas for the radial Dirac equation.) An- 
other method is to apply the general theory. It proves [2J that for the case of 
rational Weyl functions the corresponding operators with difference kernels 
can be inverted explicitly by formulas from [I]. 

The case of the skew-self- adjoint Dirac type system with the rational 
Weyl function was treated in [21J . It was shown that any strictly proper 
rational p x p matrix function is the Weyl function of a skew-self-adjoint 
Dirac type system on semi-axis and the solution of the inverse problem was 
constructed explicitly similar to the self-adjoint case treated in [20J. 

The analogues of the operators with difference kernel for the skew-self- 
adjoint system (II. ip are bounded operators Si in Lp(0, 1) (0 < / < oo), which 
have the form [301 E3] 

Sif = Sf = f(x) + llf +t t k { r + X 2 ~ t ) k { r + l~ X ) m drf(t)dt, (1.3) 

where sup 0<a:< ^ ||A;(a;)|| < oo. The kernel of the operator Si is denoted by K: 

1 f x+t , (T + x — t\ , (T + t — x\ * , 

K(x ' t:) = M,_ t| *(— 2— M— 2— ) *• (L4) 
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In this paper we show that for a Weyl function of the form 

<p(X) = <f>(X) exp{-2iXD}R, D>0, (1.5) 

where is a strictly proper rational p x p matrix function, D is a p x p 
diagonal matrix, and R is a p x p unitary matrix, the corresponding operator 
S is semiseparable. Using results on the inversion of the semiseparable oper- 
ators, the inverse problem to recover the system from (p is solved explicitly. 
Putting D = 0, we get the subcase of rational Weyl functions. Some defini- 
tions and results for the general type (non-explicit) case of inverse problem 
including Theorem 12.31 and the important formula (13. 9p are also new. The 
semiseparable matrices and operators are actively studied (see, for instance, 
[TTl [T9l [HI HQ]), and their application to inverse problems is of interest, too. 

Various definitions and results on a general type inverse problem for the 
skew-self-adjoint Dirac type system and on explicit solutions of the inverse 
problem, when the Weyl functions are proper rational, are given in Section 
EJ Some properties of the operator Si of the form (11.31) are studied in Section 
[3j The explicit solution of the inverse problem for the Weyl functions of the 
form (jl.5p is contained in Section |H 

We denote by R the real axis, by R + the positive semi-axis, by C the 
complex plane, and by C + the open upper halfplane > 0. The class of 
bounded linear operators acting from Hi into H 2 is denoted by {Hi, H 2 }, 
the identity operators are denoted by /, and spectrum is denoted by o. 

2 Inverse problem. Preliminaries 

First, normalize the fundamental solution u(x, A) of system (11.11) by the initial 
condition 

u(0,X) = I 2p . (2.1) 

If 

sup ||v(x)|| < M, (2.2) 

0<x<oo 

the unique pxp Weyl matrix function <p(X) of the skew-self-adjoint Dirac type 
system (11.11) on the semi-axis [0, oo) can be defined (see also [51 |2"T| I3"3"]) 
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by the inequality 

POO 

/ [ v (\y I p ] u(x, X)*u(x, A) 
Jo 

which holds for all A in the halfplane < — M < 0. Under condition (12.21) 
such a Weyl function always exists. 

Consider the case of the so called pseudo-exponential potentials |21j, 
which are denoted by the acronym PE. A potential v GPE is determined by 
three parameter matrices, that is, by the n x n matrix a (n > 0) and two 
p x n matrices 6\ and 9 2 , which satisfy the identity 

a -a* = i{9 1 9* l + 9 2 9* 2 ). (2.4) 

The pseudo-exponential potential has the form 

v{x) = 29* 1 e ixa *i:(x)- 1 e ixa 6 2 , (2.5) 

where 

E(x) = I n + [ A(t)jA(t)*dt, A(x) = [ e' ixa 9 1 e ixa 9 2 ] . (2.6) 
Jo 

By Proposition 1.4 in [21] the pseudoexponential potential v, i.e., the poten- 
tial given by (12. 5p is bounded on the semi-axis. The Weyl function of the 
system (11 .11) with v GPE is a rational matrix function, which is also expressed 
in terms of the parameter matrices |21j : 

<f(X) = i9* 1 (XI n - /3)- 1 9 2 , (3 = a- i9 2 9* 2 . (2.7) 

In spite of the requirement 

P-F = i(e 1 e* 1 -e 2 e^, (2.8) 

which is implied by the equalities (12.41) and f3 = a — i9 2 9 2 , any strictly proper 
rational matrix function can be presented in the form (12. 7p . The inverse 
problem to recover v from the strictly proper rational matrix function (p is 
solved explicitly in [21], using a minimal realization of <p and formula (I2.5p . 



tp(\) 



dx < oo, 



(2.3) 
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When ( 12.21) is true, inequality (12. 3p implies other inequalities: 



sup 

x<l,Q\<-M 



e txA u(x, A) 



¥>(A) 



< oo for all < / < oo, (2.9) 



which can be treated as a more general definition of the Weyl function. 

Definition 2.1 Let the system be given on the semi-axis [0, oo). Then 
a p x p matrix function y?(A) analytic in some halfplane SA < — M < is 
called a Weyl function of this system, if inequalities Ii2.9\) hold. 



If 



sup ||u(x)|| < oo for all < Z < oo, (2.10) 

0<x<l 



then there is at most one Weyl function. 

Definition 2.2 The inverse spectral problem (ISP) for system U.l\) on the 
semi-axis is the problem to recover v(x) satisfying \2.9\) and li2.1U\) from the 
Weyl function if. 

For an analytic matrix function ip satisfying the condition 

sup ||A 2 (</?(A) -a/A) || < oo, (2.11) 

GA<— M 

where a is some p x p matrix, the solution of the inverse problem always 
exists (see Lemma 1 [30J for the scalar version of this result and the matrix 
case can be proved quite similar). 

The general (non-explicit) procedure to solve ISP is described in [2SJ [2H1 
|3(H [33] . Fix a positive value I (0 < / < oo). The first step to solve ISP is to 
recover a p x p matrix function s(x) with the entries from L 2 (0, 1) (I < oo), 
i.e., s(x) e Lp Xp (0,Z) via the Fourier transform. That is, we put 

s(x) = ^e-^U.m.^ [ e^'A-V(A/2K (A = £ + irj, r? < -2M), 

27r J-a 

(2.12) 

the limit l.i.m. being the limit in L 2 (0, 1). As (12.121) has sense for any I < oo 
the matrix function s(x) is defined on the non-negative real semi-axis x > 0. 
Moreover, it is easily checked that s is absolutely continuous, it does not 
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depend on the choice of rj < —2M, s' is bounded on any finite interval, and 
s(0) = 0. To define the operator Si we substitute k(x) = s'(x) into (11.31) . 
Next, denote the p x 2p block rows of u by U\ and lo 2 \ 

uj 1 (x) = [I p 0}u(x,0), u 2 (x) = [0 l p ]u{x,0). (2.13) 

It follows from ([HI]) that u{x,0)*u{x,0) = I 2p . Hence, by flUD and fl2~T3l) 

we have 

v(x)=w' 1 (x)u 2 (x)*, (2.14) 
and uji, u 2 satisfy the equalities 

ui(0) = [I p 0], u x ul = I p , ou[ul = 0, uiuZ = 0. (2.15) 

It is immediate that uj\ is uniquely recovered from oj 2 using (I2.15p . 
Finally, we obtain uj 2 via the formula 

uj 2 (1) = [0 I p ] - J (Sfh'ixj) [I p s(x)}dx (0 < I < oo), (2.16) 

where Sf 1 is applied to s' columnwise. 

From the considerations in [29, 30J (see also similar constructions in [31], 
where the Weyl theory for the linear system auxiliary to the nonlinear optics 
equation is treated) it follows that one can solve ISP under requirements on 
ip and s(x) weaker than (12. lip . Namely, we assume 

sup HA)||<oo, (2.17) 

QX<~M 

ip(X) G Lp Xp (— oo, oo), A = £ + 277 (-oo < £ < oo) for all r] < -M, (2.18) 
s(0) = 0, sup ||jfe(aO|| < oo for all < I < oo, k(x) := s'(x), (2.19) 

0<x<l 



e' cx \\k(x)\\dx < oo (2.20) 



for some c > 0. 



Theorem 2.3 Let the matrix function ip be analytic in the halfplane < 
—M and satisfy the relations (2.11) and $2.180 . Let also the matrix func- 



tion s(x) defined via ip by formula $2.120 be absolutely continuous and satisfy 
$2.19) and $2.29) . Then ISP has a unique solution, which is given by formu- 
las $2.140 - $2.16) . where Si > I has the form fOj) with k = s' . 
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3 Factorization of S and operator identity 



Consider again the operator S = Si. It is easy to see that functions, which 
are bounded on the interval, can be approximated in the L 1 -norm by the 
continuous functions. As k = s' is bounded on the finite intervals, one can 
see that the kernel K of S, which is given by (11.41) . is continuous with respect 
to x and t. Hence, the kernel of Sf 1 is continuous with respect to x, t, and 
I ([23], p. 185). Therefore, Sf x k has the form ^S^kjfa) = k(x) + ki(x), 

where ki is continuous, and the matrix function ISf^kjty) is well-defined: 



S^k^l) = k(l) + h(l) = k(l) + lim ((S^k)(x) - k{x)^j (0 < / < oo). 

(3.1) 

To express v in terms of (^Sf 1 ^ (I) we need some preparations. According 
to [31] there are triangular operators Vi G (Lp(0,Z), 1/2(0, 1)}, such that 



(V l f)(x) = f(x)+ V-(x,t)f(t)dt, ViAVr 1 = iuj^x) u^ty-dt, 
Jo Jo 

V-(x, t) does not depend on /, and the operators V\ and Vf l map functions 
with bounded derivatives into functions with bounded derivatives. Moreover, 
as bounded functions on an interval can be approximated in the L 1 -norm by 
the continuous functions, it follows from the construction in |34J that V(x, t) 
(x > t) is continuous with respect to x and t. 
Next, introduce the operator 

(V l f)(x) = f(x) + J*V.(x-t)f(t)dt, 7_W:=^(^Vi)W, (3.3) 

where u\\ is the first p x p block of uii, and put 

V r .= = 1+ [ V.(x,t) ■ dt. (3.4) 
Jo 

It is easy to see that V\A = AVi, and so the second equality in f)3.2p yields 

ViAV' 1 = iwi(x) / ut(ty ■ dt. (3.5) 
Jo 
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By (12.131) we see that u>n(0) = I p . Hence, using definition (13.31) one gets 

W p = I p + J X (jj-(Vr x wn)){x-t)dt (3.6) 

= I P + £ (^(\r 1 u> n ))(t)dt= 

Formula (I3.6P implies V l ~ 1 cu n = I p . Moreover, from [291 E3] it follows that 
under the conditions of Theorem 12.31 the equalities 

(VTV)(x) = [/p s(x)} (3.7) 

and 

sr 1 = v:vi (3.8) 

are also true. 

Remark 3.1 Under the conditions of Theorem \2.3\ we have 

v(l) = (SfV)(l). (3.9) 

Indeed, using Ii3.8\) and changing variables I and x into x and t, correspond- 
ingly, we rewrite 112.16)) in the form 

oo 2 (x) = [o i P ]-J o {yj)(tyv x [i P s(t)]dt. (3.io) 

As V- does not depend on I we have (y x s')(t) = (V/s')(t) for t < x < I. 
Thus, according to ^3. 7| ) and I13.1U\) . we get 

w' 2 (x) = -(V l s')(x)*u 1 (x). (3.11) 

Multiplying both sides of Ii3.11\) by uo\ from the right and taking into account 
{2. 1J$ and Ii2.15\) . one derives —v(x)* = — (VJs')(x)* ; i.e., the equality 

v(x) = {V lS ')(x) (3.12) 

is true. As V(x, t) is continuous, taking into account \3. 3\) and ft3.4\ ) we see 
that (Vis') (x) — s'(x) is continuous. It is also immediate from Ij3.4\ ) that 

(VCf)(x) = f(x) + [ l V^t,x)*f(t)dt. (3.13) 



Hence, according to A3. A3. 8\) . and A3. ISA) we have 

(sr l s')(l) = (v lS i)(l). (3.14) 

Finally, formula ( TO|) follows from A3. ISA) and j3.14\) . 
By (ESK the equality 

as - sa* = vr 1 (wr 1 - {v^vf 1 )*) {vr 1 )* 

is valid for S = Si. Therefore, taking into account (13.51) and (13.71) one can 
see that 5* satisfies the operator identity 

AS-SA* = iUU*, n=[$! $ 2 ], ^g = g, <5> 2 g = s(x)g. (3.15) 

Here G {C p , Lp(0, /)} {k = 1,2) and C denotes the complex plane. This 
identity differs from the identity AS - SA* = i(^^* 2 + $ 2 ®1) [351 US] for 
an operator with difference kernel. Matrices satisfying a discrete analogue of 
(I3.15P were treated in [17] . The operator identity (I3.15P for the case, when 
k in (11.31) is a vector, was studied in |25j. It could be useful also to prove 
(I3.15P directly. In fact, we prove below a somewhat more general identity. 



Proposition 3.2 Let the operator S in L"i(0,l) (0 < I < oo) be defined by 

sf = f{x) + ^/^K 111 ^)^! 111 ^)^^' (3 - 16) 

where sup 0<x< i (||/c(x)|| + ||A;(a;)||) < oo. ThenS satisfies the operator identity 

AS-SA* = i[ (I p + if)(x)f(t)) ■ dt, (3.17) 
Jo 

where ip(x) = J* k(t)dt, ip(x) = J* k(t)dt. 

Proof. Using (13.161) and changing the order of integration we have 

ASf = Af + i f 7l {x,t)f{t)dt, (3.18) 
J o 

7l(M ) : = i r r\r-±j^y k r-±ifJL\ driy . (3 . 19 ) 



2 Jo J\y-t\ ^ - 



Taking into account that for the scalar product (•,•)/ in Lp(0,l) we have 
(A* f, g)i=(f, Ag)i, rewrite SA* in the form 

SA*f = A*f-i [ 72 (x, t)f(t)dt, (3.20) 
Jo 
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First, consider the case t > x. From (13.191) . after changes of variables £ = 
(r + y — t)/2 and r\ = t — y + £, we get 

7l(M) = i r r k (i+iz±y k (i±pi) drdv ,3.22) 



2 Jo Jt-y V 2 / V 2 

/ [ V k(0Kt - y + O^dy = [ [ k(Z)k(t - y + Z)dyd£ 
Jo Jo Jo Js, 

[ k(0Hv)dvdt 

o Jt-x+i 

Next, calculate 72(2, t) (t > x). From (13.211) it follows that 

7 2 (x,t) = 721 (x, t) +722 (a;, t), (3.23) 

where 

m ~\[ C^'^^P^**'' (3 - 24) 

*(«,*) ~\££kp^)kp^)M*. (3.25) 

Replace the variable r by 77 = (r + y — x)/2 in (I3.24p , then change the order 
of integration, and after that put £ = x — y + 77 and change the order of 
integration again to obtain 

72iOM)= / f* k(Qk(ri)drid£. (3.26) 
Jo Jo 

In (13.251) . replace r by £ = (r + x — y)/2, change the order of integration and 
put r) = y — x + £. We get 

l22(x,t)= / / k(0k(v)dvd£. (3.27) 



J § 
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By (133211 . (EE23J , dSSHD, and (1X27|) the equality 



7i(x,t)+ 7 2(x,t) = / k{i)k{r 1 )dr ] di = ij{x)^{t) (3.28) 

is true for t > x. Using similar calculations one can show that ( 13.281) holds 
also for x > t, i.e., ( I3.28P is true for all < x, t < I. Finally, formulas (I3.18p . 
( KM . and (I3T28D yield (IHTTjl . ■ 

4 ISP and semiseparable operators 5/ 

In this section we consider matrix functions of the form 

<p{\) = z0*(AJ n - P)- 1 9 2 e- 2iXD R, (4.1) 
D = diag{<ii, . . . , d p }, d kl > d k2 > for &i > fc 2 , (4.2) 

where 9j (j = 1, 2) is an n x p matrix with the m-th column denoted by 6j )7n , 
f3 is an n x n matrix, R is a p x p matrix, and D is a p x p diagonal matrix. 
We do not suppose here that 9j and (3 satisfy the identity (12. 8p . 

Proposition 4.1 Let matrix function tp be given by jjj4.1\ ). Then, the matrix 
function s, which is defined via (p by ( \2.12\l . has the form s = CR, where 
C = [ci C2 . . . c p ] , the columns c m (p > m > 1) being given by the formulas 

c m {x) =0 for < x < d m , (4.3) 

rx~d m 

c m (x) = 261 / exp{2it(3}dt6 2) m for x > d m , (4.4) 
Jo 

and the function ip is the Weyl function of system U.l\) with potential v 
satisfying h2.10\) . 

Proof. First, choose M > such that a(/3 + iMI n ) C C+, where a means 
spectrum and C + is the open upper halfplane. According to (14. ip <p(\) is 
analytic and the function A(/?(A) is bounded in the halfplane < —M. 
So, the conditions (12.171) and (I2.18P on <p are fulfilled. The fact that s is 
absolutely continuous and satisfies conditions (I2.19P and (12.201) is immediate 
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from (14.31) and (14.4)) . Therefore, after we have proved (14. 3)) and (14.41) , it will 
follow from Theorem 12.31 that ip is the Weyl function of system (11.11) with 
potential v satisfying (I2.10p . 

Now, let us prove (14.3)) and (14.41) . As \(f(\) is bounded, one can rewrite 
point wise limit: 

s=[ Cl c 2 ... c p ]R, c m {x) = —6{ / e **(*-*») (4.5) 
x\-\\I n -2(3)- 1 d£6 2 , m (\ = Z + i v , V <-2M). 

Introduce the counterclockwise oriented contours, where f may take complex 
values: 

r+ = [~a, a] (J{£ : If I = a, 3£ > 0}, r~ = [-a, a] \J{£ : |f | = a, 3£ < 0}. 

For A = f + 277 and for the fixed values of 77 < — 2M, it follows from (14.51) 
that 



Cm. \ X I 



■-9{ lim / e iX( - x - dm) \- 1 (\I n -2(3)- 1 d£6 2 , m (4.6) 



in the case x > d m , and 



C77?. ( ^£ ) 



-0* lim / e iX{x - dm) \~ l (\I n - 2(3)- l d£0 2 , m (4.7) 

7T a^oo J r - 



in the case x < d m . As e* A ^ _dm ^A _1 (A/ ri — 2/?) -1 is analytic with respect to 
f inside T~ and on the contour itself, equality (14. 3)) is immediate from (14. 7p . 

Next, consider the > d m . For sufficiently large a all the poles of 

(XI n — 2(5) ~ l (and the pole f = — it] of A -1 ) are contained inside r+ and 
taking into account (14 .6p we have 



Cm. I J 



■-61 f e^ x - d ^\- 1 (XI n -2(3r l d^ m . (4.8) 



Let us approximate /? by matrices f3 e such that — /3 e || < e and det (3 £ 7^ 
(if det /3 7^ we put (3 = (3 £ ) . It is easy to see that 

\-\Mn - 2(5 £ )- 1 = (2/5 £ )- 1 ((A/ n - 2&)- 1 - X-'ln). (4.9) 
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For sufficiently small e all the poles of (XI n — 2/? e ) 1 are contained inside r+ 
and we have 

/ e**(AI„ - 2/3 £ )~ 1 d£ = e"* exp(2ix/3 £ ) (x > 0). (4.10) 
Finally, using (O)) - (O0j) we get 

rx~d m 

c m {x) = \im.29\ I exp{2it(3 £ }dt8 2m . 
Jo 



Hence, formula (I4.4p is immediate. ■ 

Remark 4.2 Note that the matrix functions tp of the form ( [^. 1\) in general 
position do not satisfy A2.11\) and so they do not satisfy in a scalar case 
conditions of Lemma 1 l3(Jf . but the conditions of Theorem \2.3\ are fulfilled. 

By Proposition 14.11 the matrix function k in the expression (11.41) for the 
kernel of the operator Si, generated by the Weyl function ip of the form ( 14. ip . 
is given by the formula 

k( x ) = s'(x) = 2B\e 2 ^v X [x)R, v : = {exp(-2id m (3)e 2im } p m=1 , (4.11) 

X(x) = dmg{xi(x),X2(x),...,Xk(x)}, X m(x) := I °' ^ x < d m, 

According to ( 11 Ah and (14.111) we have 

K(x,t) = 261 J ex P (i(r + x - t)(3^jQ(r, x,t) exp - i(r + t - x)(3*^jdr6 1 , 

(4.12) 

where 

Q(r,x,t) = vx(- ± ^)RR*x(- ± ^y- (4.13) 

The matrix function Q(r,x,t) is piecewise constant with respect to r and 
without loss of generality we assume Q(0, x, t) = 0. It is easy to see that 
Q(r, x, t) has only a finite number of jumps {Qj}. Moreover, if <r(/3) n<j(/3*) = 
0, the matrix identity i(/3Xj — Xj/3*) = Qj always has the solution Xj. 
Therefore we have 

e irP Q,e- irP * = ^(e^Xje-^y (4.14) 
13 



Hence, according to (14.121) and (14. 14j) we can express the kernel K (x, t) of S 
explicitly in terms of matrix exponents and {Xj}. It follows also from (11.41) 
that 

K{x,t) =K{t,x)\ (4.15) 

and so we need to simplify (14.121) only for x > t. 

Another approach to the presentation of K in terms of matrix exponents 
is given in the following lemma. 

Lemma 4.3 Put 

9j (r):=[0 e^V*' 
Then we have 

(^)W=e^Q,e-^. (4.17) 
Proof. By (I4.16P we have 

^-gj = iflgj + [0 e^]E 3 e^ 

and (I4.17P is immediate. ■ 

Recall [18] that the operator S is called semiseparable, when K admits rep- 
resentation 

K(x,t) = F x (x) Gx(t) forx>t, K(x, t) = F 2 {x)G 2 {t) for x < t, 

(4.18) 

where F\ and F 2 are px p matrix functions and G\ and G 2 are p x p matrix 
functions. For the operator S to be semiseparable, assume 

RR* = I p . (4.19) 

Then the matrix function Q has the form 

Q(r,x,t) = ^x( r + 2~ X )^ for x>t ( 42 °) 



Ip 






, Ej : = 








-i(3* 
Qj —ip 



(4.16) 



= iPg j + e irP Q j e-*P -ip 9i , 



14 



Rewrite (jOjl as 

D = diag{di/ pi , . . . , d k I Pk }, p 1 + . . . + p k = p, d kl > d k2 > for k x > k 2 , 

(4.21) 

and put 

Q 3 = isPjv*, P j = diag{0, . . . , 0, I Pj , 0, . . . , 0}. (4.22) 



Remark 4.4 Some notations. Further we consider K(x, t) (x > t) on the 
intervals d m < t < min(x, d m +i), where we choose such m that the in- 
equalities d m < x hold. If d\ > we put d = and include the interval 
do < t < min(x,di) into consideration. If x > d k , we include the interval 
d k < t < d k +i (d k +i = x). Some matrix functions, like B{t) and C{t), will 
be considered on the intervals as above, but with x = I. In the following, 
in all such cases we simply write d m < t < d m+1 . We also assume that 
Y^jLi • • • = 0, when m = 0. 

Proposition 4.5 Let the matrix function tp be given by where D sat- 

isfies Iji4.21\ ) and R is unitary. Assume also that the matrix identities 



i((3X J -X j (3*)=Q 3 , (4.23) 

where Qj are given by {^.22$ , have solutions Xj . Then the operator S , which 
is defined via tp by formulas < fi.3|) . k = s' and Ii2.12\) . is semis eparable, and its 
kernel K(x, t) (0 < x, t < I) is given by relation ( [^.15] ) and by the equalities 

K{x, t) = 261 (e 2ixf3 Z m e- 2UI3 * - e 2i ^Z m ^ 9 1 (d m < t < d m+1 ) (4.24) 

for t < x < I. Here 

m m 

Z m = Y,Xj, Z m = J2( e ^P( 2i d^))x j e^p(-2id j (3*). (4.25) 

3=1 3=1 



Moreover, there are self-adjoint solutions of {4-23 ) and we suppose Xj = X* 
in UJSD and jJJty . 
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Proof. First, note that we can choose Xj = X* because the adjoint of each 
solution of (I4.23P also satisfies (14.231) . 

Next, using (Q2j) . (14^201) . and ffl~22j) we get equalities 



" b px+t 

K(x, t) = 291 { ^2 / _(expi(r + x-t)p)Qj 

^ j =1 Jx-t+2dj 

x exp ( — i(r + t — x)(3*)dr\0i 



(4.26) 



for t < x and d m < t < d m+x . From fl4TT4|) and ffl~26D it follows that (1421 
holds. Formula (14.151) was derived earlier. ■ 



Remark 4.6 By \2.lJ$ - $2.l6t) and \J^.l 1\ ) the equality v(x) = is valid for 
< x < d\ in the case d\ > 0. This fact corresponds to the inequality 



sup 





e e J 


- ^x)e~ 2aD ' 




M 









< oo, 



(4.27) 



which can be easily checked directly and is implied also by Ii2.9\) . 

When the operator S = I + J ' K(x, t) ■ dt is semiseparable and its kernel K 
is given by f)4.18p . the kernel of the operator T = S~ l is expressed in terms 
of the 2p x 2p solution U of the differential equation 



G 2 {x) 



(— U) (x) = H(x)U(x), x>0, U(0) = I 2 p. 

where 

H(x) := B(x)C(x), B(x) -- 

Namely, we have (see, for instance, [T8] ) 

T = S- 1 = 1+ [ T(x,t) ■ dt, 
Jo 

f C(x)U(x)(l 2p ~-P*)U(t)-'B(t), x>t, 
[X, > \ -C(x)U(x)P x U(t)- 1 B(t), x<t. 



(4.28) 



, C(x) = [ F 1 {x) F 2 (x) ] . 

(4.29) 



(4.30) 
(4.31) 
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Here P x is given in terms of the p x p blocks Uzi(l) an d ^22(0 °f 







(4.32) 



and the invertibility of is a necessary and sufficient condition for the 

invertibility of S. 

If K admits the representation 



K(x,t) 



Ce xA (l 2 p- P)e- tA B, x>t, 
-Ce xA Pe~ tA B, x<t, 



(4.33) 



where A, B, and C are constant matrices, then U is calculated explicitly 
In our case a representation 



K(x,t) 



C m t xA {l2 v — Pm) e ^Bm, t < X < I, d m < t < d m+ i, 
C m & Pm& B mj X *C t *C Z, G? m <C X <C (i m _)_i, 



(4.34) 



A = 2i 



(4.35) 



where p = n and 

" /3 
/3* 

easily follows from (14.151) and (14.241) . However, (14.341) is insufficient for the 
explicit construction of U and we shall construct U and T explicitly, using 
more general formulas (I4.28l) - (l4.32p . For this purpose we introduce B(x) and 
C(x) (0 < x < I) by the equalities 



B(x) = V2 



-2*0 Z 



-2ixf3* 



-2ix/3* 



m < x < d m j r \ / 



(4.36) 



C(x) = V261 [ e 2i *P e 2 ^Z m - Z m e 2 ^* ] (d m < x < d m+ i), (4.37) 
where Z m = and Z m = Z^ are defined in (14.251) . 

Proposition 4.7 Let the conditions of Proposition \4^5\ be fulfilled and let S 
be defined via by formulas $1.3\) . k = s' and $2.12\) . Then the operator 
T = S' 1 is given by formulas $JJU\) - $J?%\) , fflTZFfy , and 



U(x) = n m e- xA e xA ^U{d m ) (d m <x< d m+1 ), U(0) = I 2n , (4.38) 
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where A is defined by {4-3 5 ) and 
A x ■= A + 2Y Y 



^1 , 



j _ 7 



7 
In 



0i9l 







/'/I 



Moreover, we have 

U(x)*JU(x) = J, [/(a;)" 1 = JU(x)*J*, J :-- 







(4.39) 



(4.40) 



(4.41) 



Proof. Recall that B and C are recovered from K by formulas (14.181) and 
(T4T29]) . In view of (1415]) and (Oil we have 



F 1 (x) = V20{e 2lxP , F 2 (x) = Gi(a?) 



(4.42) 



G^z) = v^^e-^'-e- 24 ^)^ (d m < a; < d m+1 ), (4.43) 
G 2 {x) = Ft(x)*. (4.44) 



Therefore, formulas (14.291) and (I4.42I) - (I4.44I) imply that B and C correspond- 
ing to S are given by (1436]) and (14371) . It follows from (14^29]) and (jQ5j) - 
f l437j) that 



ff(ar) = 2n m e- xA Y m e xA n^ (d m < x < d m+1 ] 
where Y m is given in (I4.39j) . Q m is given in (14.401) . and 

nz 1 -- 



T 7 



o /„ . 

According to (TOED, (j4\35]) . and (Q5]) we get 



(4.45) 



(4.46) 



d 

fix 



[/ ) fa:) = Q m e~ 



xAl A X 



•A™ .A ) e 



[/(d m ) = H{x)U{x) 



for d m < x < d m+ i, and so £/ of the form (14.381) satisfies (14.281) . In other 
words, formulas (I4.36l) - (l4.38p define explicitly B, C and U, which are used 
in the expressions (I4.3ip and (I4.32p to construct the kernel of T = S~ l . 
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It remains to prove (14.411) . Note that 

J A* J* = -A, JQ* rn J* = ft" 1 , JY^J* 
Hence, we have JH*J* = -H, i.e., 

4-(U(x)*JU(x)) =0. 
dx 



-Y 



Formula (I4.4ip follows from (I4.48P and from U(0) = I 2n - 
Taking into account (I4.36I) - (I4.38I) we get 



F{x) := C(x)U(x) = V29l[I n - Z m ]e x ' 



G(t) := Uit^Bit) = v / 2?7(rf m )- 1 H m e- t - 4 " 



U(d m ), 
0i, 



dm < x < d 



m+l) u r7i 



d m < t < d. 



(4.47) 



(4.48) 



(4.49) 
(4.50) 



m+l • 



Corollary 4.8 Let the conditions of Proposition \4-5\ be fulfilled. Then the 
kernel T(x, t) of the operator T = Sf 1 has the form 



T(x,t) 



F(x)(l 2n -^P*)G{t), x>t, 
-F(x)P x G(t), x<t, 



(4.51) 



where F and G are given by fi4.49j ) an d ft4.50\ ). 

By Km . fl439|) . and KM) for d m < t < d m+1 we get 



G{t)k(t) = V2U(d rn )- 1 E rn e- tA -(2Y m )e tA 
/— ../T , i_ d 





' In ' 








m 

uJ2 p j R - ( 4 - 52 ) 



From Remark 13. II and formulas (14.301) . (14.511) . and (I4.52p the explicit solution 
of the inverse problem is immediate. 
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Theorem 4.9 Let the Weyl matrix function ip be given by ( |^.ip 7 where D 
satisfies and R is unitary. Assume also that the matrix identities 

^4-23^ , where Qj are given by \-j.22\j , have solutions Xj = X* . Then the ISP 



solution v is given by the formula 

V{1) = k{l) + F{l){hn ~ P X ) 



N 



(4.53) 



x e 



m=l 



) 


" In ' 








3=1 



where k is given by {J^.ll ), U is given by \4.38 ), P x is given by l[4-32\ ), and 



S m is given by lji4-40 )- The number N in the sum is chosen in the following 
way: if I < d\ then N = 0; if dj < I < dj + i then N = j ; if I > dk then 
N = k. We put d m = d m for m < N and dx+i = I- 
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